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We  examined  die  spin-up  from  rest  of  water  in  a  rectangular  cyclinder.  The  presence 
of  corners  in  the  cylinder  causes  the  formation  of  eddies.  We  found  that  the  number  of  ed¬ 
dies,  as  well  as  eddy  size,  position,  and  rotation  rate  were  dependent  on  the  aspect  ratio  of 
the  cylinder,  the  depth  of  the  fluid,  and  the  Anal  angular  velocity  of  the  cylinder.  Two  time 
scales  were  found  to  be  important  in  this  experiment:  the  traditional  Ekman  number  based 
on  depth ,  which  defines  the  time  scale  required  for  spin-up  and  an  additional  Ekman  number 
based  on  the  cylinder  length  which  provides  some  information  about  the  evolution  of  the 
fluid  pathlines  in  route  to  spin-up.  This  second  Eckman  number  appears  to  provide  an  ex¬ 
planation  for  both  the  agreement  and  disagreement  of  the  experimental  results  herein  and 
previously  published  results. 
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Chapter  1 

INTRODUCTION 

Spin-up  is  the  flow  of  a  contained  fluid  that  changes  from  one  state  of  rigid  body 
rotation  to  another  when  the  container’s  rotation  rate  is  increased  from  a  steady,  initial  angu¬ 
lar  speed  of  to  a  steady,  final  angular  speed  of  Qf.  Applications  of  spin-up  include  mo¬ 
tions  in  the  Earth’s  core,  the  stability  of  liquid-filled  projectiles,  the  transient  development 
of  wind-driven  coastal  upwellings,  and  the  interior  rotation  of  the  sun  (Benton  and  Clark, 
1974).  Spin-up  is  also  encountered  in  models  of  global  ocean  circulation  that  are  started 
from  rest  and  attempt  to  develop  a  quasi-steady  circulation  pattern  (Pond  &  Pickard,  1983, 
pi 83- 190).  These  problems  all  involve  nonlinear  systems  whose  multi-parameters  are  dif¬ 
ficult  to  identify  and  isolate.  Herein,  we  consider  a  seemingly  simple  system  that  is  in  fact 
somewhat  complicated  -  the  spin-up  from  rest  of  water  in  a  right,  rectangular  cylinder.  The 
flow  depends  on  a  number  of  parameters  including  water  depth,  H,  horizontal  aspect  ratio, 
6,  and  Qf.  We  describe  an  experimental  investigation  of  the  effects  of  varying  these  parame¬ 
ters  on  the  evolution  of  the  flow  during  spin-up  and  the  pathlines  of  the  flow  at  spin-up. 

Despite  the  complications  that  are  specific  to  any  given  system  undergoing  spin-up, 
some  of  the  basic  physics  is  similar.  In  particular,  the  flow  involves  viscous  boundary  layers 
at  the  walls  and  an  inviscid  flow  in  the  interior.  The  prototype  system  for  study  is  the  spin-up 
of  water  in  a  right  circular  cylinder.  As  the  angular  velocity  of  the  cylinder  is  increased  from 
Qo  to  Qf,  the  interior  fluid  is  unaffected  with  pressure  gradient  forces  balancing  centrifugal 
forces.  However,  because  of  the  no-slip  boundary  condition,  the  fluid  on  the  boundaries  has 
the  angular  speed  of  the  cylinder.  The  resulting  large  viscous  stresses  rapidly  cause  the  fluid 
near  the  cylinder  boundaries  to  accelerate.  In  the  absence  of  horizontal  boundaries,  the  angu- 
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lar  velocity  at  the  sidewall  diffuses  into  the  interior  due  to  viscosity.  The  presence  of  horizon¬ 
tal  boundaries,  though,  creates  a  different  phenomenon.  As  the  angular  velocity  of  the  fluid 
in  the  horizontal  boundary  layers  increases,  the  centrifugal  force  soon  overcomes  the  pre¬ 
vailing  pressure  gradient  force,  causing  the  fluid  to  be  thrown  radially  outward.  Conserva¬ 
tion  of  mass  requires  that  the  fluid  thrown  outward  be  replaced  by  fluid  from  the  inviscid 
interior  and  a  much  slower,  secondary  flow  is  established.  This  secondary  flow  is  the  key 
to  the  spin-up  process.  The  secondary  flow  convects  high  angular  momentum  fluid  away 
from  the  cylinder  sidewalls  toward  the  rotation  axis.  As  the  fluid  moves  into  the  interior, 
its  angular  momentum  is  conserved.  Thus  the  secondary  flow  transports  higher  angular-mo¬ 
mentum  fluid  into  the  cylinder  interior  faster  than  would  occur  if  viscous  diffusion  was  the 
only  process  involved.  Once  the  fluid  is  rotating  uniformly  at  the  new  angular  velocity,  Qf, 
we  say  it  is  spun-up. 

The  physics  described  above  applies  to  spin-up  in  the  rectangular  cylinder  in  a  gener¬ 
al  sense;  however,  the  presence  of  corners  causes  flow  separation  that  results  in  the  formation 
of  eddies  that  evolve  within  the  interior  of  the  flow.  Numerous  and  varied  eddy  configura¬ 
tions  result  from  the  spinning-up  of  fluid  in  a  rectangular  cylinder.  Let  n  be  the  integer  value 
nearest  to  the  value  of  the  horizontal  aspect  ratio,  6.  In  general,  we  have  found  that  if  n  is 
odd,  then  the  pathlines  of  the  flow  at  spin  up  form  6  circular  cells,  which  are  antisymmetric 
about  the  diagonals  of  the  cylinder  and  have  vorticity  (before  spin-up)  at  the  narrow  walls 
that  has  the  sense  of  the  cylinder’s  rotation.  This  result  is  in  contrast  with  that  of  van  Heijst, 
Davies,  and  Davis  (1990),  hereafter  vHDD,  who  observed  that  the  vorticity  (before  spin-up) 
at  the  narrow  walls  has  the  opposite  sense  of  the  cylinder’s  rotation.  If  n  is  even,  the  flow 
evolves  into  an  odd  number  of  cyclonic  and  anti-cyclonic  cells  whose  centers  are  not  on  a 
common  horizontal  axis.  This  result  is  true  when  the  change  in  rotation  rate,  AQ,  is  relatively 
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small  and  the  water  is  relatively  deep.  However,  as  we  increase  AQ,  or  decrease  the  water 
depth,  H,  the  flow  patterns  change.  For  example,  our  results  mimic  those  of  vHDD  in  experi¬ 
ments  for  which  an  Ekman  number  based  on  the  (long)  length  of  the  cylinder  is  small  enough, 
hi  Chapter  5,  we  describe  the  results  of  changing  6,  H,  and  Qf. 

A  more  thorough  review  of  previous  research  concerning  spin-up  will  be  presented 
in  Chapter  2.  Chapter  3  describes  the  basic  equations  governing  the  flows  and  a  linear  solu¬ 
tion  for  spin-up  in  a  circular  cylinder.  The  experimental  equipment  and  procedures  are  dis¬ 
cussed  in  Chapter  4,  and  the  results  are  presented  in  Chapter  5. 
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Chapter  2 

Literature  Review 

Previous  research  on  spin-up  has  consisted  primarily  of  theoretical  and  experimental 
studies  of  spin-up  in  axisymmetric  cylinders.  We  briefly  review  this  work  to  discuss  the  ba¬ 
sic  physics  of  fluid  undergoing  spin-up.  (See  Greenspan  (1968)  and  Benton  &  Clark  ( 1974) 
for  comprehensive  reviews .)  We  then  discuss  previous  work  on  spin-up  in  non-axisymmet- 
ric  cylinders. 

The  first  linear  model  describing  the  spin-up  of  a  homogeneous  fluid  in  a  closed 
axisymmetric  container  was  provided  by  Greenspan  and  Howard  (1963).  They  examined 
the  limiting  case  of  an  infinite  circular  cylinder  of  radius  A.  In  this  case,  they  found  that  the 
mechanism  for  spin-up  is  the  diffusion  of  vorticity  that  arises  at  the  lateral  boundaries.  Rigid 
body  rotation  is  achieved  in  a  dimensional  time  of  order  Tv  =  A2/v,  where  v  is  the  kinematic 
viscosity  of  the  fluid.  Next,  they  analyzed  the  opposite  extreme,  that  of  two  parallel  coaxial 
infinite  disks  separated  by  a  distance  of  2H.  In  this  case,  they  found  that  the  spin-up  con¬ 
sisted  “of  three  distinct  phases:  the  development  of  the  Ekman  layer,  the  inviscid  fluid  spin- 
up,  and  the  viscous  decay  of  residual  oscillations.”  They  determined  that  the  dimensional 
spin-up  time,  which  Greenspan  and  Howard  defined  as  the  time  required  for  the  interior  flow 
to  approach  the  final  state  of  rigid  rotation  within  e“l,  is  of  order  Te=  (H2/v£2)1/2.  Using 
these  results  for  two  parallel  coaxial  infinite  disks,  they  used  a  perturbation  expansion  in 
half-powers  of  E  to  study  the  case  of  an  arbitrary  axisymmetric  closed  container.  Greenspan 
and  Howard  noted  that  if  the  container  has  a  vertical  sidewall  a  viscous  boundary  layer,  with 
a  double  structure  different  from  that  Ekman  boundary  layer,  is  formed.  These  layers  affect 
only  a  small  portion  of  the  fluid  near  the  sidewall  prior  to  spin-up.  Thus,  the  spin-up  time 
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scale  is  still  of  order  Te  unless  the  container  is  elongated  in  the  vertical  direction  to  the  extent 
that  the  ratio  H/L  (depth/width)  s  E-1^4  in  which  case  the  time  scale  for  spin-up  is  again 
the  viscous  time-scale,  Tv.  Greenspan  and  Howard  also  discuss  the  two  effects  of  having 
a  free  surface.  First,  the  Ekman  layer  at  the  top  boundary  is  removed  so  the  secondary  flow 
changes  from  a4-cell  circulation  to  a  2-cell  circulation  (See  Figure  2.1).  Second,  they  note 
that  the  pressure  gradient  that  arises  from  the  parabolic  shape  of  the  free  surface  induces  a 
radial  motion  toward  the  interior  in  addition  to  the  secondary  flow  produced  by  the  Ekman 
layer  on  the  bottom  of  the  cylinder.  This  radial  motion  acts  to  decrease  the  spin-up  time. 
The  Froude  number,  F  =  Q2L2/gH,  where  £2L  is  the  characteristic  velocity,  is  a  measure  of 
the  relative  importance  of  this  ffee-surface  effect. 
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Figure  2.1  Cross-sectional  View  of  the  Secondary  Flow,  (a)  rigid  lid  (b)  free  surface. 


The  first  non-linear  model  of  the  spin-up  of  a  homogeneous  fluid  within  a  circular 
cylinder  of  finite  length  and  radius  was  done  by  Wedemeyer  (1964).  He  divided  the  flow 
into  two  distinct  regions,  the  boundary  layer  flow  near  the  end  walls  of  the  cylinder  and  the 
core  flow.  He  derived  a  nonlinear,  inhomogeneous  hyperbolic  partial  differential  equation 
for  the  radial  component  of  the  core  flow.  The  solution  to  the  equation  showed  that  the  core 
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flow  is  divided  into  two  regions:  an  inner  region  in  which  the  fluid  particles  do  not  rotate, 
and  an  outer  region  in  which  the  particles  spin-up.  The  regions  are  separated  by  a  cylindrical 
front  of  fluid  which  moves  toward  the  axis  of  rotation. 

Weidman  (1976)  did  extensive  experimental  work  to  examine  the  theoretical  work 
of  Wedemeyer  (1964).  Weidman  noted  that  one  result  of  the  Wedemeyer  theory  is  an  0(1) 
velocity  discontinuity  that  arises  in  the  inviscid  Wedemeyer  model  when  the  full  nonlinear 
Ekman  suction  is  included  in  the  spin-up  equations.  Weidman ’s  measurements  of  the  azi¬ 
muthal  velocity  for  spin-up  in  a  circular  cylinder  indicated  that  this  velocity  discontinuity 
is  smoothed  out  in  a  shear  layer  whose  thickness  is  a  function  of  radius  and  time  and  scales 
with  E^H.  He  noted  that  the  velocity  distributions  changed  depending  on  whether  the  ratio 
of  the  wall  acceleration  period,  Ta  =  l(S2j  _  Qf)/al,  where  a  is  the  cylinder  acceleration  rate, 
to  the  ordinary  Ekman  spin-up  time,  Tg,  was  less  than  or  greater  than  0(1).  Weidman  ob¬ 
served  that  for  small  accelerations,  Tj/Ie  >  0(1),  measured  velocities  agreed  well  with  the 
Wedemeyer  theory  away  from  the  propagating  shear  layer.  The  effect  of  the  shear  layer  is 
local  and  does  not  affect  the  large  scale  flow.  At  large  accelerations,  T/Te  <  0(  1 ),  Weidman 
found  significant  deviations  between  the  experiments  and  Wedemeyer ’s  theory.  In  this  case, 
the  effect  of  the  shear  layer  was  not  negligible. 

As  noted  before,  the  models  described  above  were  all  for  axisymmetic  geometries. 
A  study  of  spin-up  phenomena  in  non-axisymmetric  cylinders  was  done  by  van  Heijst 
( 1989).  He  presented  experimental  observations  of  spin-up  for  several  geometries:  (i)  an 
annular  region  between  two  coaxial  cylinders  with  a  radial  banrier  between  the  cylinder 
walls;  (ii)  a  semicircle  region;  (iii)  a  circular  region  with  a  radial  barrier  extending  from 
the  center  to  the  tank  wall;  and  (iv)  a  region  enclosed  by  two  non-concentric  cylinders  con¬ 
nected  by  a  radial  barrier.  vHDD  (1990)  studied  spin-up  in  a  rectangular  container.  In  each 
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of  these  experiments,  the  authors  observed  three  main  stages  in  the  spin-up  process:  (i)  the 
starting  flow,  characterized  by  zero  absolute  vorticity,  (ii)  flow  separation,  and  (iii)  a  subse¬ 
quent  organization  of  the  flow  into  an  array  of  alternately  cyclonic  and  anticyclonic  circular 
cells.  In  the  latter  experiments,  vHDD  noted  that  the  central  cell  in  the  free-surface  experi¬ 
ments  was  always  cyclonic. 
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Chapter  3 

THEORETICAL  CONSIDERATIONS 

In  this  chapter,  we  present  the  equations  of  motion  associated  with  the  rectangular 
cylinder  as  well  as  the  boundary  and  initial  conditions  associated  with  the  problem  of  spin- 
up.  We  then  describe  a  solution  for  the  flow  at  some  small  time  increment,  t  =  e.  Lastly, 
we  review  the  linear  solution  for  spin-up  in  a  circular  cylinder. 

3.1  Equations  of  Motion  for  Soin-up  in  a  Reciamular  Cylinder 

A 

Consider  a  fluid  of  depth  H  in  a  rectangular  cylinder  rotating  at  -  Qk,  with  period 
T  =  2n/Qf.  In  the  following  analysis,  we  use  a  rectangular  coordinate  system  (x,  y,  z)  that 
is  rotating  with  the  cylinder.  The  origin  is  at  the  water  surface  in  the  center  of  the  cylinder 
as  shown  in  Figure  3.1. 


Figure  3.1  Schematic  View  of  the  Rectangular  Cylinder 

The  Navier-Stokes  equations  for  a  flow  with  velocity  vector  u *  =  (u*,  v*.  w*>  are  (the  *  de¬ 
notes  dimensional  values): 


V* • u  =  0 


3.1 
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fijj£  -  2 at  x  „•  =  -  ivy  -  **  +  vV'V  3.2 

The  boundary  conditions  are  u*  =  0  at  x*  =  ±  ^L;  v*  =  0  at  y*  =  ±  jW;  and  w*  =  0  at  z*  = 
0,  and  z*  =  -H  (where,  for  simplicity,  we  consider  the  free  surface  to  behave  as  a  rigid  lid). 
The  pressure,  p*.  consists  of  perturbations  p*  from  the  hydrostatic  pressure. 

P*  =  -  Q*8*z*  +  P  ( x,y,z,t )  3.3 


To  non-dimensionalize  the  above  equations,  we  scale  the  x,  y,  and  z  directions  by  their 
lengths,  L,  W,  and  H  respectively,  and  time  by  the  final  rotation  rate  Qf  Since  we  are  inter¬ 
ested  in  an  unsteady  flow  due  to  a  change  in  rotation  speed.  A,  we  scale  velocities  by  the 
appropriate  length  scale  times  A .  Thus,  we  have: 


3.4 


t  = 


P  g£2fH2A 


3.5 


Then  the  appropriate  non-dimensional  numbers  are  the  Ekman  number,  E  h  v/QfH2,  the 
Rossby  number,  e  =  A/Qf,  and  aspect  ratios,  a,  0,  and  8  defined  as  (H/L)2,  (H/W)2,  and  L/W 
respectively.  So  the  incompressible  Navier-Stokes  equations  in  dimensionless,  coordinate 
form  are: 


du.dv.dw_r. 

to  +  Ty+te-° 


3.6 


dll  -/„du  i  „dM  ,  dw  .  2v  _  dp  r,  d2u  ,  od2u  ,  d2u  ,  ,  - 
~dt  +  e{udx  +  vdy  +  W~^}  ~~  ~aj;  +  E{a—  +P—  +  3-7a 


dzJ 


dx 


dx2  H  dy2  dz2  ‘ 


dv  ,  _/„dv  ,  „dv  ,  „,dv  .  _  adp  c/„d2v  ,  od2v  .  d2v ,  , 

Ji  +  cluai  +  ''Ty  +  wTzl-2iu-  ~^  +  Elal?  +  ^  +  ^>  37b 
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£  +  +  w&l  =  -  f  +  El&$  +  +  &Z)  3.7c 

dt  dx  dy  dz  dz  dx2  dy2  dz2 

The  boundary  conditions  ucu  =  v  =  w  =  0  at  x  =  —1;  y  ~  —A ;  and  z  =  -1. 


3.1  Streamlines  for  the  Initial  Flow 


If  we  define  a  stream  function  as 

u  =  —  ipy  and  v  =  \j)x  3.9 

then  the  vorticity,  <v,  is  V2^.  At  the  onset  of  spin-up,  the  inner  flow  is  quiescent  relative 

to  the  inertial  reference  frame,  while  the  flow  on  the  boundaries  rotates  with  the  boundaries. 
In  the  rotating  frame,  the  sign  of  the  vorticity  in  the  interim*  is  opposite  that  of  the  cylinder’s 
rotation  (i.e.,  anticylonic)  and  is  zero  on  the  boundaries;  therefore,  Vfy  =  2  with  tp  =  0  on 
the  boundaries.  Hence,  a  boundary  value  problem  at  the  onset  of  spin-up  is  (following 
vHDD) 


VV  =  2  =  3l0 

y  =  0  at  y  =  -ii 

A  particular  solution  of  (3.10)  is  xpp  =  (y  +  1/2)2  -  (y  +  1 12)  =  y2  -  1!4\  then  (3.10)  can  be 
transformed  into  the  following  homogeneous  form  of  Laplace’s  equation: 


v2n  =  0  V>  =  \/4-y2atx=  3U 

V-Oaty-  -i± 

The  solution  to  (3.1 1)  is  found  using  separation  of  variables  and  fourier  series  to  be: 


g  A  sin(2n  +  l>r(y  +  ;)  sinh(2/i  +  l)^(j  -  x)6  +  sinh(2n  +  l)rf.x  + 
jr3^  (2  n  +  l)3  sinh(2n  +  l)xd 


-w-1 


3.12 
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A  plot  of  tp  and  its  contours  is  shown  in  Figure  3.1.  vHDD  observed  that  experimental 
streamlines  agreed  well  with  the  theoretical  results  described  above  for  t/T  =  .08  after  the 
onset  of  spin-up. 

3.3  Linear  Solution  for  Spin-uo  in  a  Circular  Cylinder 


It  is  useful  to  review  the  linear  solution  for  spin-up  in  a  circular  cylinder  provided 
by  Greenspan  and  Howard  (1963)  in  order  to  understand  the  mechanism  for  spin-up  and  the 
differences  between  spin-up  in  a  circular  cylinder  and  spin-up  in  a  rectangular  cylinder. 
Consider  a  fluid  in  a  circular  cylinder  of  radius,  A,  and  depth  H,  rotating  in  rigid  body  rota¬ 
tion  about  its  center  with  constant  angular  velocity  £2j.  Consider  a  cylindrical  coordinate  sys¬ 
tem  ( r,  0,  z)  rotating  with  the  cylinder  whose  origin  is  in  the  center  of  the  cylinder.  For  sim¬ 
plicity,  we  consider  the  fluid  to  have  a  rigid  lid.  We  seek  a  solution  for  the  azimuthal  velocity 
in  the  interior  of  the  cylinder. 


The  dimensional  Navier-Stokes  equations  for  this  axially  symmetric  flow  with  ve¬ 
locity  vector,  ii*  =  (u*.  v*,  w*),  and  perturbation  pressure,  p ,  are: 


1  d(rV)  |  dw*  _  0 


*  -x  * 

r  dr 


dz * 


3.13 


Oil  _  )il  -  2Qv*  =  - 


-  o  +  V(A*U*  -  4^) 
Dt  r  Q  dr*  r  2 


3.14a 


Dv*  ,  u*v* 


Dt' 


T  +  ~~  +  2 Qu=v{A*v*  -  ^ ) 


3.14b 


Dw* 

Dt* 


{dp  *  . 

--—  +  vA„ 


3.14c 


If  we  suppose  the  angular  velocity  is  impulsively  increased  by  a  small  amount,  e£2,,  then 

*  *  * 
the  boundary  and  initial  conditions  are  u  =  eQk  x  r  on  the  solid  boundaries,  and  u  = 
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Ofort  s  0.  To  non-dimensionalize  the  above  equations,  we  introduce  the  following  non-di¬ 
mensional  variables 


=  ^z  =  - 
H,  H' 

A* 

P=  ' 


t  =  t*Q, 


ii  =  _«_ 

HeQ' 


and  E  — 


3.15 


3.16 


QeQH2'  QH2' 

where  £  =  (Qf  -  Qi)/Q;  is  the  Rossby  number  and  E  is  the  Ekman  number.  Considering  the 

case  of  £  <3  1,  we  can  neglect  the  nonlinear  terms.  The  linear  initial-boundary  value  problem 

becomes 


1  d(ru)  ,dw  =  n 
r  dr  d: 


%  -  1  +  ECJu  -  4) 


dp 


dt 


dr 


dv 


^  +  2«  =  EGdv  -  -4) 
dt  rl 


3.17 


3.14a 


3.14b 


dw 

dt 


dp 

-  -f  +  EAw 
dr 


3.14c 


with  initial  condition  u  =  0  for  t  <  0  and  boundary  condition  u  =  k  x  r  on  all  solid  bound¬ 
aries. 


Define  a  stream  function,  ip,  such  that 

u  =  -  jxpz,  andw  =  \xpr  319 

where  the  subscripts  indicate  partial  differentiation.  The  coordinate  equations  for  (3.18)  in 
the  case  of  an  axisymmetric  cylinder  become 

-  jxpa  -  2v  =  -  pr  +  E(A{-  j\pz)  +  ^y),  3.20a 

vt  -  jV>z  =  +  E{Av  - 

r* 


3.20b 
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and  y\pn  =  -  Pz  +  EA{j\pr) .  3.20c 

Rearranging  (3.20b)  and  eliminating  p  from  (3.20a  and  c)  provides 

£(Jv  ~  pj)  -  v,  +  jipz  =  0,  and  3.21a 

E{A(—  ftyzz  iy^Pr)r)  3-  -^3*)  ~f^Prr  ~  f^Pr  3"  ^Pzz)t  2vz  —  0  .  3.21b 

rhh 

The  initial  condition  is  tp  =  v  =  0  for  t  <■  0,  and  the  boundary  condition  is  xp  =  ^  -  0, 
and  v  =  r  on  the  lateral  boundary. 

Since  the  time  scale  of  interest  for  the  infinite  parallel  disks  was  shown  to  be  order 
E~1/2Q  we  begin  our  analysis  of  the  circular  cylinder  by  scaling  as  follows: 

x  =  E\  and  xp  =  Efy,  3.22 

where  we  have  scaled  the  radial  and  vertical  flow  to  be  higher  order  in  Ekman  number  than 
the  azimuthal  flow  which  is  0(1).  Substituting  r  and  x  into  equations  (3.21)  we  obtain 

£5(dv  -  ^)  -  vr  +  |xr  =  0,  and  3.23a 

jXzz  -  (7Xr)r)  +  *f)  +  f  fcrr  “  jXr  +  Xzz)z  +  2vz  =  0  .  3.23b 


The  flow  can  be  divided  into  an  interior  and  boundary  layer  flow  such  that 

V  =  V7  +  V8,  and*  =  xl  +  XB  •  3.24 

The  equations  (3.23)  for  the  interior  flow  to  0(E°)  become 

4  ~  jXz  =  0,  and 


4  =  0. 


3.25a 

3.25b 


Thus, 


15 


v7  =  vV.rJ.and 

3.26a 

Z7  *  f  Vr  +  • 

3.26b 

This  defines  the  interior  flow  in  terms  of  (as  yet)  undetermined  functions,  v7  and  Xq,  of  r 
andr. 

In  the  boundary  layer,  we  introduce  a  stretched  vertical  coordinate,  such  that 

z  =  -  ^  +  in  the  lower  boundary  layer,  and  3.27a 

2  -  2~  in  upper  boundary  layer  .  3.27b 


Then  to  0(E°)  the  boundary  layer  equations  are 


Lower 

+  7*f  -  0, 

and  tX^k  -  2v®  =  0; 

3.28a 

Upper 

'ft  -  7*f  =  ft 

and  jXBm  +  2vf  =  0; 

3.28b 

The  boundary  conditions  for  t,  -*  oc  are  v8  —  xB  =  0;  for  £  =  0,  they  are  v5  +  v7  =  r, 
Xfl+X/  =  0,and^  =  O(£1/2). 

Taking  first  integrals  of  (3.28a),  we  have 

vf  +  jXB  =  /for),  and  jjgg  -  2V5  =  g(r,r);  3.29 

where  /and  g  are  equivalent  to  constants  of  integration.  Eliminating  xB  from  (3.29)  gives 

vBm  +  4vB  =  0  .  3.30 

The  boundary  conditions,  vB(oo)  =  0,  v*(0)  =  r  -  v1,  and  v^l0  =  0,  imply  that 


v5  =  (r  -  vfy,  r))e  ^  cos  ^  . 


3.31 
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A  similar  analysis  shows  that  the  same  formula  applies  on  the  upper  boundary.  From  (3.28) 


XB(0)  =  ^r(r  -  v7)  on  z  =  -  and  3  32a 

*B(0)  =  -  ^r(r  -  v7)  on  z  =  ± .  3.32b 

The  boundary  condition,  XB  +  X1  =  0.  (3.32),  and  (3.26b)  imply 

^v£(r,r)  -  ^(r.r)  =  ~  v7),  and  3.33a 

Jv7(r,  r)  +  x{>(r,r)  =  ^  (r  -  uf0) .  3.33b 

Adding  (3.33a  and  b),  we  get 

vi(r,r)  =  2(r  -  v7) .  3.34 

Thus,  the  interior  azimuthal  flow  is 

vV.r)  =  r(l  -  e~2x) .  3.35 

°r  vV.fl  =  rd-e-vy  336 


From  this  result,  we  see  that  the  interior  azimuthal  flow  is  linearly  dependent  on  r  and  the 
spin-up  time  scale  is  of  order  Tg  =  (H2/v&)  • 
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Chapter  4 

Experimental  Equipment  and  Procedures 

The  experiments  were  conducted  in  the  Department  of  Mathematics  Fluid  Mechan¬ 
ics  laboratory.  The  laboratory  equipment  consists  of  a  real-time  computer,  the  turntable 
and  drive,  two  flow  visualization  systems,  and  the  rectangular  tank  with  prepared  water.  Fig¬ 
ure  4.1  is  a  schematic  of  this  system,  which  we  describe  below. 


Figure  4.1  Schematic  of  the  Computer,  Turntable,  and  Drive  Assembly 

4.1  Real-Time  Computer 

The  command  signal  to  the  motor  drive  was  from  a  Digital  Equipment  Corporation 
Micro  VAX  Workstation  II  with  analog-to-digital,  and  digital-to-analog  capabilities,  as 
well  as  two  real-time  clocks.  In  order  to  minimize  transients,  we  wrote  a  program  to  ramp 
the  voltage  smoothly  from  a  user-specified  initial  voltage  to  the  desired  final  voltage.  The 
program  allows  the  user  to  input  the  slope  of  the  ramp  (i.e.,  the  time  to  reach  the  final  volt- 
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age).  For  the  experiments  described  herein,  we  used  a  1  sec  ramp,  sending  1000  incremental 
voltages  to  the  controller  at  1000Hz.  The  final  voltage  was  then  sent  to  the  controller  at 
SOOHz  for  as  long  as  desired.  To  determine  the  relationship  between  voltage  and  the  rotation 
rate,  Q,  of  the  table,  we  measured  the  time  for  5  rotations  of  die  table  at  12  different  voltages. 
The  calibration  curve  is  shown  in  Figure  4.2.  The  curve  is  approximately  linear  on  a  log-log 
scale.  Our  operating  window  is  between  0  and  8  volts,  corresponding  to  rotation  rates,  Q, 
from  0  to  2.1  rad/s. 


OQ 


O 


3.5 


3 


2.5 


“°-5 I 0.5  1  •  1.5 


1.5 


log  (voltage) 


Figure  4.2  Motor  Calibration  Curve 


42  Turntable  and  Drive 

The  rotating  table  is  a  Lintech  model  MR-308180.  The  platform  is  8  inches  in 
diameter  and  is  mounted  to  the  table  shaft  with  4  screws  located  11/16  inches  from  the  center 
of  the  table  (see  Figure  4.3).  A  large  square  piece  of  white  lucite,  13  1/3  inches  on  a  side, 
is  mounted  on  the  platform  using  2  screws  located  3  inches  from  the  center  of  the  table.  The 


Tank 


Pine 


Table 


Figure  4.3  Schematic  of  the  Table/Tank  Assembly 
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rectangular  cylinder  is  set  on  the  table  using  a  pin  and  bracket  system,  which  centers  the  cyl- 
inder  to  within  1/500”.  The  pin  has  a  plug  on  one  end  that  sits  in  the  center  of  the  table; 
the  other  end  fits  tightly  into  a  hole  on  the  bottom  of  the  tank.  Four  metal  brackets  hold  the 
container  in  place.  Measurements  with  a  Brown  and  Taylor  dial  indicator  show  that  the  table 
itself  rotates  about  its  center  with  an  accuracy  of  1/500”.  The  table  is  level  to  within  5.8  x 
10"6  radians. 

The  cylinder  exhibits  a  high-frequency  vibration  that  varies  in  frequency  and  inten¬ 
sity  with  table  speed.  This  vibration  causes  high-frequency  waves  on  the  water  surface. 
Figure  4.4  shows  the  frequency  spectra  of  these  waves  versus  table  speed.  The  vibration  is 
minimal  when  the  system  turns  without  the  platform,  lucite  plate,  or  tank  mounted.  It  inten¬ 
sifies  as  each  of  these  components  is  mounted;  apparently  the  vibration  increases  as  weight 
is  added  away  from  the  table  center.  We  placed  a  rectangular  piece  of  lead  between  the  bot¬ 
tom  of  the  turntable  and  the  laboratory  bench  in  an  attempt  to  minimize  this  vibration.  We 
also  isolated  the  entire  rotating  table  assembly  from  extraneous  laboratory  vibrations  by 
mounting  it  to  a  structural  wall  of  the  laboratory. 

The  table  drive  is  provided  by  a  DigiPlan  Electronic  Motion  Control  BL  Brushless 
Servo  Drive  System.  The  system  consists  of  a  BL150  Drive  with  velocity  feedback,  a  TO170 
transformer,  and  an  ML-3450,  dc  brushless  motor.  The  signal  travels  from  the  computer  to 
the  drive,  through  the  transformer  where  it  is  converted  from  120V  AC  to  170V  AC,  back 
to  the  drive,  and  to  the  motor.  The  velocity  feedback  signal  then  comes  back  from  the  motor 
to  the  drive  which  corrects  the  velocity  to  ensure  steady  motion.  This  servo  system  broke 
down  when  the  input  voltage  was  greater  than  about  8  V  (2. 1  rad/s)  at  which  point  the  rota¬ 
tion  rate  of  the  motor  was  unsteady.  We  note  that  prior  to  the  acquisition  of  this  motor  we 
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(d)  0.52  rad/s,  (e)  0.65  rad/s,  (f)  0.78  rad/s,  (g)  1.04  rad/s,  (h)  1.30  rad/s,  (i)  1.56  rad/s,  and  (j)  1-82  rad/s.  The  ordinates  ore  ln(a). 
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drove  the  rotating  table  with  a  stepper  motor.  Qualitatively,  all  of  the  results  we  report  herein 
were  present  when  we  used  the  stepper  motor,  including  the  high-frequency  vibration. 

43  Flow  Visualization 

We  visualized  pathlines  of  the  fluid  using  both  a  digital  imaging  and  an  analog  video 
system.  The  first  is  a  Kodak  EktaPro  EM  Motion  Analyzer,  Model  1012.  This  system  con¬ 
sists  of  a  high-gain  imager  with  a  Nikon  35mm  lens,  a  video  processor,  and  a  Sony  Trinitron 
Color  Video  Monitor.  The  imager  sensor  has  a  192  X  239  pixel  array  with  256  grey  shades. 
The  processor  can  store  1637  full  frames  of  digital  images  in  electronic  memory  at  various 
rates  from  50  to  1000  Hz.  To  obtain  lower  recording  rates,  and  correspondingly  longer  time 
series,  one  must  trigger  the  imager  externally.  We  used  a  Tectronix  TM502 A  signal  genera¬ 
tor  to  provide  a  0-5V  square  pulse  at  frequencies  between  15-50  Hz  for  this  purpose.  The 
processor  has  a  reference  reticule  with  built-in  X-Y  electronic  crosshairs  and  a  reference 
marker.  We  mounted  the  imager  on  a  Bogen  3051  professional  tripod.  The  imager  views 
the  rotating  tank  through  a  standard  mirror  mounted  at  about  a  45°  angle  above  the  tank. 

We  used  the  EktaPro  to  track  the  location  of  the  eddy  centers  as  a  function  of  time. 
(See  §5 .5)  To  this  end,  we  used  the  reticule  to  determine  the  position  of  the  center  of  the  eddy. 
The  errors  inherent  in  this  procedure  are  due  to  (1)  the  subjectivity  of  the  observer,  (2)  the 
limited  spatial  resolution  of  the  imager,  and  (3)  the  optical  set-up  of  the  imager/mirror  sys¬ 
tem.  By  comparing  the  results  obtained  from  two  observers,  we  found  that  the  errors  due 
to  ( 1)  are  small  compared  to  those  due  to  (2)  and  (3).  The  errors  due  to  (2),  which  is  the  dis¬ 
tance  between  pixels,  introduce  a  standard  deviation  in  position  location  of  ±.07  in.  To  deter¬ 
mine  the  errors  due  to  (3),  we  tracked  a  mark  on  the  tank  bottom  during  rotation.  Figure 
4.5  displays  the  pathline  we  obtained.  The  error,  which  is  a  combination  of  error  due  to  (2) 


24 


and  (3),  is  ±.23  in.  We  transferred  the  images  to  VHS  tapes  using  a  Panasonic  AG- 1960 
Multiplex  VCR  linked  to  the  processor. 

The  analog  video  system  was  a  Sony  Hi-8  Video  Camera  Recorder,  Model 
EVO-150TR.  It  records  at  30  Hz  directly  onto  tape.  We  mounted  the  camera  directly  above 
the  rotating  tank  and  used  a  remote  controller  for  direct  recording  of  the  flow  paths.  We 
played  the  tapes  on  a  Sony  EVO-9800  VCR.  To  obtain  “instantaneous”  images  of  the  flow, 
we  photographed  the  monitor  using  a  Nikon  FM2  camera  with  a  50mm  lens.  These 
photographs  are  shown  in  §5.1  -5.5.  We  found  two  problems  associated  with  photographing 
the  monitor.  The  first  is  that  the  monitor  resolution  limits  the  photograph  resolution.  The 
second  is  that  the  curvature  of  the  monitor  distorts  the  image.  This  distortion  is  especially 
apparent  in  Figures  5 .5.  It  introduces  errors  of  ±0.20  inches  into  the  measurements  discussed 
in  §5.6. 

4.4  The  Rectangular  Tank  and  Water 

We  constructed  the  1 1  x  3|  in2  tank  out  of  clear  lucite,  fastening  the  with  a  clear, 

fast-curing  solvent  cement.  To  vary  the  aspect  ratio  (length/width)  from  three  to  two,  we 
used  quarter-  and  half-inch  inserts  that  fit  tightly  inside  the  box.  The  tank  was  filled  with 
distilled  water  to  the  desired  depth.  The  depth  was  measured  using  a  Lory  Type  C  point 
gauge  accurate  to  0.01  mm.  Fluid  motions  were  observed  using  red,  green,  and  blue  food 
coloring.  A  drop  of  the  food  coloring  was  applied  to  the  surface  and  generally  sank,  creating 
a  column  through  the  water  and  a  blob  on  the  bottom.  In  most  of  the  figures  in  §5, 3  drops 
were  applied  symmetrically  along  the  two  long  sidewalls. 
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Figure  4.5  Data  Points  Showing  Measurement  Error  for  the  Position  of  the  Eddy  Center 
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Chapter  5 
RESULTS 

In  this  chapter  we  present  and  discuss  our  observations  of  the  pathlines  of  fluid  that 
spins  up  in  a  rectangular  cylinder  when  the  rotation  rate  is  impulsively  changed  from  Q,  to 
Qf.  We  begin  by  discussing  a  specific  “comparison”  case  for  which  the  horizontal  aspect 
ratio  of  the  cylinder,  6  =  3,  Q0  =  0,  Qf  =  0.52  rad/s,  and  the  water  depth,  H  =  5  cm.  This 
experiment  is  similar  to  one  by  vHDD,  to  which  we  compare  our  results.  Next,  we  discuss 
qualitatively  the  effects  of  changing  6,  Qf ,  and  H.  Lastly,  we  show  measurements  of  the 
path  of  an  eddy  center  and  a  separation  point  in  the  flow.  In  all  of  these  cases,  we  assess 
the  flow  as  an  observer  in  the  inertial  reference  frame  (not  in  the  rotating  frame)  looking 
down  on  the  tank.  (We  did  several  experiments  in  the  rotating  frame  and  our  conclusions 
did  not  change.)  Figure  5. 1  is  a  schematic  of  the  cylinder  plan  view;  it  shows  what  we  will 
refer  to  as  the  “leading”  comers,  A  and  C,  the  “trailing”  comers,  B  and  D,  the  distances 
referred  to  as  the  “length”  and  “width”  of  the  tank,  and  the  angle  0  that  we  refer  to  in  §5.2. 


Figure  5.1  Plan  Mew  of  the  Rectangular  Cylinder 
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5.1  Comparison  Case 

The  series  of  photographs  in  Figure  5.2  show  the  pathlines  of  the  flow  for  an  experi¬ 
ment  in  which  the  fluid’s  angular  velocity  was  changed  from  zero  (initially  it  is  at  rest)  to 
Qf  =  0.52  rad/s  (the  final  state  is  rigid  body  rotation)  when  6  =  3  and  H  =  5cm.  We  chose 
this  experiment  as  a  baseline  for  comparison  with  our  other  experiments  and  with  an  experi¬ 
ment  performed  by  vHDD.  In  the  following,  we  describe  horizontal  flow  observations,  ver¬ 
tical  flow  observations,  and  compare  our  results  to  the  flow  patterns  of  a  similar  experiment 
performed  by  vHDD. 

5.1a  Horizontal  Flow  Observations 

When  the  tank  began  to  spin,  the  fluid  in  the  center  of  the  tank  was  stationary.  Fig 
5.2a  shows  that  within  three  seconds,  a  large  anticyclonic  flow  developed  that  was  similar 
to  the  initial  flow  described  in  §3.2  in  the  rotating  reference  frame.  The  fluid  moved  quickly 
along  the  x-sidewalls  from  D  to  A  and  B  to  C.  In  accordance  with  the  no-slip  boundary 
condition,  the  fluid  on  the  wall  stayed  on  the  wall,  so  the  edge  of  the  anticyclonic  cell  pro¬ 
vided  some  indication  of  the  boundary  layer  thickness.  Within  5-10  seconds,  four  cyclonic 
eddies  formed  in  the  comers  of  the  tank.  Figure  5.2b  shows  that  the  two  cells  in  the  leading 
corners  were  noticeably  larger  than  those  in  the  trailing  comers.  By  about  1  1/3  rotations 
(Figure  5.2c)  the  leading  comer  cells  grew  while  the  trailing  corner  cells  were  no  longer  vis¬ 
ible.  Figure  5.2c  also  shows  a  separation  point  in  the  flow,  with  two  small  counter-rotating 
eddies  on  either  side.  This  separation  point  remained  a  distinct  feature  of  the  flow;  we  show 
measurements  of  its  path  in  §5.6.  As  the  cyclonic  cells  grew,  their  centers  moved  out  of  the 
leading  comers.  We  show  measurements  of  such  a  path  in  §5.5.  Figure  5.2d  shows  that  the 
cyclonic  cells  divided  the  initial  anticyclonic  flow  into  a  primary  anticyclonic  cell  in  the 
center  of  the  tank,  and  two  secondary  anticyclonic  cells  in  the  trailing  comers  of  the  tank. 


t&siSat'  h->- 1 
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Figure  5.2e  shows  that  by  2  rotations  of  the  cylinder,  the  cyclonic  cells  spanned  the  width 
of  the  tank  and  the  center  anticyclonic  cell  was  more  nearly  circular.  A  sharp  edge  was  appar¬ 
ent  between  the  two  cyclonic  and  the  center  anticyclonic  cells.  The  distinction  was  not  so 
clear  at  the  y-sidewalls  of  the  box,  where  the  secondary  anticyclonic  cells  were  still  present. 
As  the  flow  continued  to  spin  up,  the  cyclonic  cells  continued  to  grow  with  their  centers  mov¬ 
ing  just  past  the  y-centerline  of  the  tank;  the  secondary  anticyclonic  cells  occupied  the  entire 
y-sidewall  with  the  prominent  portion  of  the  eddy  now  in  the  leading  corners  of  the  tank. 
Figure  5.2f  shows  the  spun-up  fluid.  The  pathlines  formed  an  essentially  3-cell  pattern, 
with  the  center  anticyclonic  cell  flanked  by  two  cyclonic  cells.  The  center  of  the  anticyclonic 
cell  formed  on  the  axis  of  rotation  and  remained  there  throughout  spin-up.  The  size  and 
position  of  the  secondary  anticyclonic  cells  after  spin-up  depended  on  6,  Qf,  and  H  as  will 
be  discussed  in  §5.2  -  5.5.  The  flow  remained  roughly  antisymmetric  about  a  horizontal  di¬ 
agonal  of  the  cylinder  throughout  the  spin-up  process. 

5.1b  Vertical  Flow  Observations 

A  definitive  determination  of  the  vertical  flow  was  difficult  to  obtain  for  the  experi¬ 
ments  herein.  Since  the  dye  is  heavier  than  the  water,  a  slow  sinking  of  dye  was  indistin¬ 
guishable  from  a  slow,  secondary,  downward  flow.  However,  a  secondary  upward  flow  was 
unmistakable.  From  repeated  observations,  an  approximate  graph  of  the  vertical  flow  in  the 
primary  eddies  was  obtained  (See  Figure  5.2.1).  Figure  5.2.  la  shows  that  the  dyelines  in  the 
primary  anticyclonic  eddy  spiraled  up  and  into  the  interior  of  the  eddy.  The  vertical  coordi¬ 
nate  increased  like  ea0  where  a  is  0(0. 1),  while  the  radial  component  decreased  like  Re~P0, 
where  R  is  the  cell  radius  and  j3  is  0(0.01).  Vertical  motion  in  the  center  of  the  anticyclonic 
eddy  appeared  to  be  downward .  Figure  5.2.1b  shows  that  dyelines  in  the  two  cyclonic  eddies 
spiraled  up  and  away  from  the  center  of  the  eddy.  The  vertical  coordinate  increased  like  ea0, 


Figure  5.2.1  Approximate  Trace  for  the  Observed  Vertical  Flow  Patterns 
the  Eddies,  (a)  anticyclonic,  (b)  cyclonic. 
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where  a  is  0(0. 1),  while  the  radial  component  increased  like  Re^0  -  R,  where  R  is  die  eddy 
radius  and  0  is  0(0.1).  The  vertical  motion  in  the  center  of  the  cyclonic  cells  was  upward. 

5.1c  Comparison  with  vHDD’s  Experiment 

The  flow  described  above  differs  fundamentally  from  the  flow  pattern  of  a  similar 
experiment  performed  by  vHDD.  Table  1  shows  the  relevant  parameters  and  non-dimen¬ 
sional  numbers  for  the  two  experiments,  which  both  had  6  =  3  and  Q0  =  0.  (We  note  that 
vHDD’s  fluid  was  either  water  or  a  dilute  salt  solution.  It  was  not  clear  to  the  reader  which 
is  the  case  for  this  comparison  experiment.  The  value  of  v  given  in  Table  5. 1  was  obtained 
from  calculations  using  given  non-dimensional  parameters.)  The  flow  in  their  experiments 
began  as  in  our  comparison  case  with  the  initial  anticyclonic  cell  and  2  cyclonic  cells  that 
grew  out  of  the  leading  comers.  vHDD,  however,  observed  that  the  2  cyclonic  cells  then 
divided  the  original  anticyclonic  cell  into  two  anticyclonic  cells  which  remained  at  the  ends 
of  the  tank.  The  two  cyclonic  cells  in  the  vHDD  experiments  moved  toward  the  axis  of  rota¬ 
tion  in  the  center  of  the  tank  where  they  merged  into  one  cell.  The  pathlines  to  the  final  flow 
(rigid  body  rotation)  consisted  of  three  cells,  a  central  cyclonic  cell  flanked  by  two  anticy¬ 
clonic  cells.  This  final  result  is  in  complete  contrast  to  our  result  of  a  central  anticyclonic 
cell,  flanked  by  two  cyclonic  cells. 

We  did  observe  an  experiment  in  which  the  final  result  was  similar  to  vHDD’s;  al¬ 
though,  we  cannot  be  sure  that  the  evolution  to  this  state  was  the  same.  In  this  experiment, 
which  is  discussed  in  detail  in  §5 .4,  Qf  =  1 .82  rad/s.  The  relevant  parameters  for  this  experi¬ 
ment  are  also  shown  in  Table  5.1. 
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Table  5.1  Characteristic  Parameters  and  Non-dimensional  Numbers  for  One  of  v HDD’s 
Experiments  and  Two  Experiments  Herein.  The  sense  of  the  vorticity,  w,  in  the  3  main  cells 
formed  by  the  pathlines  of  the  flow  is  indicated  as  cyclonic  (+)  or  anticyclonic  (-). 


Parameters 

vHDD 

Present  Experiments 

Comparison  Case 

Large  Qf  Case 

sgn(co) 

-  +  - 

+  -  + 

-  +  - 

Qf 

0.57  s’1 

0.52  s'1 

1.82  s-1 

v 

0.005  cm2/s 

.01  cm2/s 

.01  cm2/s 

L/H 

5.9 

5.6 

5.6 

E-l/2 

160 

36 

65 

El-1'2 

949 

140 

377 

F 

.1751 

.0431 

.5277 

The  most  notable  difference  between  the  paramaters  for  VHDD’s  experiment  and 
that  of  the  comparison  case  is  in  the  difference  between  the  Ekman  numbers,  E,  of  the  two 
experiments.  The  magnitude  of  the  Ekman  number  is  affected  by  two  main  factors:  H  and 
Qf .  However,  increases  in  E-1'2  due  to  increases  in  H  did  not  affect  the  flow  characteristics 
herein.  For  example,  we  conducted  an  experiment  with  5  =  3,  Qq  =  0,  Qf = 0. 52  rad/s,  and 
H  =  10  cm.  This  experiment,  for  which  E-1/2  =  72,  showed  the  same  flow  characteristics 
as  the  comparison  case.  Increases  in  E'1^2  due  to  increases  in  Qf,  however,  did  affect  flow 
characteristics.  With  each  increase  in  Qf,  the  2  cyclonic  vortices  moved  further  toward  the 
center  of  the  tank  (See  §5.4and5.5).  Eventually,  when  Qf  =  1.82  rad/s,  the  cyclonic  vortices 
joined  and  remained  in  the  center  of  the  tank.  Based  on  this  result,  it  appears  that  the  -  +  - 
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vorticity  pattern  occurred  when  the  cyclonic  ceils  had  enough  energy  to  translate  close 
enough  to  the  center  of  the  tank  that  they  joined.  We  conjectured  that  if  the  box  length  was 
decreased,  and  the  rotation  rate  is  held  constant,  then  the  angular  momentum  of  the  comer 
cells  would  also  decrease,  and  the  cells  may  not  have  enough  energy  to  translate  to  the  center. 
Thus,  it  is  likely  that  there  are  two  scales  important  in  this  problem:  the  traditional  Ekman 
number,  E,  based  on  depth,  which  defmes  the  time  scale  required  for  spin-up  (see  §2)  and 
an  additional  Ekman  number,  El  =  v/(QfL2),  based  on  the  box  length,  which  provides  some 
information  on  the  evolution  of  the  pathlines  in  route  to  spin-up.  To  test  this  idea  further, 
we  conducted  an  experiment  in  which  £2f=  1 .82  rad/s,  6  =  3,  and  L/H  =  3.6  (the  same  parame¬ 
ters  which  gave  us  a  -  +  -  vorticity  pattern),  but  L  was  decreased  to  8  inches.  Here,  El  = 
274,  a  value  more  comparable  to  the  value  of  El  for  the  comparison  case.  As  anticipated, 
the  dyelines  in  this  experiment  returned  to  the  +  -  +  vorticity  pattern  of  the  comparison  case. 
Thus,  shortening  the  length  of  the  cylinder  decreased  the  tendency  of  the  cyclonic  cells  to 
move  into  the  center  of  the  cylinder.  Because  the  values  of  El  for  all  of  the  experiments  here¬ 
in  are  actually  about  the  same  order  of  magnitude,  further  experiments  are  needed  to  deter¬ 
mine  more  precisely  if  this  scaling  by  cylinder  length  is  an  appropriate  one. 

Another  possible  explanation  for  the  discrepancy  between  our  experiments  and  those 
of  vHDD  is  the  difference  in  Froude  number,  F,  between  the  experiments.  Greenspan  and 
Howard  ( 1963)  noted  that  the  Froude  number  is  a  measure  of  the  importance  of  free-surface 
effects  relative  to  that  of  the  Ekman  layer.  This  difference  does  not  seem  to  be  a  likely  ex¬ 
planation  as  we  conducted  experiments  with  £2f  =  1.04  and  1.30  rad/s,  discussed  in  §5.4, 
having  Froude  numbers  of  15.93  and  26.92,  respectively,  which  did  not  display  the  same 
final  result  of  a  central  cyclonic  cell  as  vHDD’s  experiment  did. 
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Another  possible  explanation  is  the  presence  of  high-frequency  waves  on  the  surface 
of  the  water  in  the  present  experiments.  vHDD  attributed  the  translation  of  the  cyclonic  vor¬ 
tex  toward  the  rotation  axis  to  the  parabolic  shape  of  the  undisturbed  free  upper  surface  of 
the  fluid.  It  is  possible  that  the  high-frequency  waves  in  the  present  experiments  act  as  a 
perturbation  of  the  free -surface  that  inhibits  the  cyclonic  vortices  from  moving  toward  the 
center  of  the  tank.  To  explore  this  possibility,  consider  the  frequency  spectra  of  the  perturba¬ 
tive  surface  waves.  The  spectral  data  for  the  Qf=  1.82  rad/s  shows  a  significant  wave  com¬ 
ponent  at  144  Hz  with  only  a  small  wave  component  at  19  Hz.  This  contrasts  with  the  data 
for  the  comparison  case,  Qf  =  0.52  rad/s,  which  had  a  number  of  small  components  across 
the  frequency  range.  This  difference  in  the  spectra  may  support  the  conclusion  that  the  sur¬ 
face  waves  affect  the  translation  of  the  vortices.  However,  the  frequency  spectra  for  Qf  = 
1.30  and  1.56  rad/s  is  similar  to  the  spectrum  for  Qf  =  1.82  rad/s;  yet,  the  dyelines  for  the 
Qf  =  1 .30  and  1 .56  rad/s  showed  the  same  central  anticyclonic  vortex  as  the  comparison  case 
while  the  Qf  =  1.82  rad/s  dyelines  showed  a  cyclonic  vortex  in  the  center. 

5.2  Effects  of  Varying  & 

A  series  of  experiments  was  conducted  to  determine  the  effect  of  decreasing  the  as¬ 
pect  ratio,  6,  from  an  approximately  odd  number,  3,  to  an  approximately  even  number,  2. 
The  first  column  of  photographs  in  Figure  5.3shows  the  dyelines  after  two  rotations  of  the 
cylinder.  Variations  in  the  aspect  ratio  did  not  significantly  affect  the  dyelines  at  this  initial 
state  of  spin-up.  The  second  column  of  photographs  in  Figure  5.3  shows  the  dyelines  after 
spin-up.  When  6  >  2.5,  the  results  were  similar  to  those  in  the  comparison  case  in  §5.1 
above.  The  anticyclonic  cells  in  the  leading  comers  of  the  cylinder  were  very  weak  and  each 
occupied  only  about  l/27th  of  the  tank.  The  centers  of  the  cyclonic  cells  moved  just  past 
the  centerline  of  the  cylinder  toward  the  trailing  comers,  leaving  the  three  main  cells  slightly 
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(a  -  g)  t/T  =  2;  (h)t/T=10,  (i)  t/T  =  11;  0)  t/T  =  10;  (k)  t/T  =  ll,(l)t/T  =  9,  (m)l/r  = 
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diagonalized.  Figure  5.3h  shows  that  when  6  =  2.38  the  secondary  anticycloruc  cells  were 
slightly  larger  and  began  to  entrain  fluid  from  both  the  cyclonic  cells  and  the  primary  anticy- 
clonic  cell.  For  6  <  25,  the  initial  state  (t fl  -  2)  was  unchanged  from  when  6  S  2.5; 
however,  the  final  state  was  different.  As  6  decreased  from  2.38  to  2,  the  two  secondary  anti- 
cyclonic  cells  became  more  prominent  and  more  off-center.  A  measure  of  the  diagonaliza- 
tion  of  the  three  main  cells  is  given  by  the  angle,  0,  between  the  line  connecting  the  centers 
of  the  three  main  eddies  and  the  y -centerline  of  the  tank.  Figure  5.4  shows  that  6  varied 
from  2°  to  30°  as  5  decreased  from  3  to  2.  When  6  =  2.24,  each  secondary  anticyclonic 
cell  was  positioned  in  the  leading  comer  of  the  tank  and  occupied  about  1/9  of  the  tank  once 
the  fluid  was  spun-up.  The  cyclonic  eddies  still  dominated,  but  0  was  increased.  When  6 
=  2.13,  each  secondary  anticyclonic  cell  grew  to  occupy  about  1/6  of  the  tank  at  spin-up. 
The  centers  of  the  anticyclonic  cells  remained  in  the  leading  comers  of  the  tank,  while  those 
of  the  cyclonic  cells  moved  toward  the  trailing  comers  of  the  tank.  Instead  of  having  an  es¬ 
sentially  3-cell  pattern  as  resulted  when  6  >  2.24,  Figures  5.3m  and  5.3n  show  that  when 
6  <  2.13  the  pattern  consisted  of  essentially  five  cells  with  a  large  anticyclonic  cell  in  the 
middle  flanked  by  a  cyclonic  and  a  anticyclonic  cell  on  each  side.  For  all  the  experiments 
above,  the  flow  remained  roughly  antisymmetric  about  a  horizontal  diagonal  of  the  cylinder. 

5.3  Effects  of  Varvine  H 

Herein  we  discuss  the  results  of  experiments  in  which  the  depth,  H,  varied  from  5 
to  0.5  cm.  The  flow  pattern  was  relatively  unperturbed  by  depth  variations  except  at  the  shal¬ 
low  depths  of  H  <1.0  cm.  In  Figure  5.5  the  first,  second,  and  third  columns  show  dyelines 
after  2, 4,  and  6  revolutions,  respectively.  The  fluid  spun-up  in  all  cases  before  6  revolutions 
of  the  tank;  when  H  =  1  and  0.5  cm,  the  fluid  spun-up  before  4  revolutions.  One  consequence 
of  decreasing  depth  was  a  corresponding  decrease  in  the  rotation  rate  of  the  cyclonic  eddies. 


(a) 


Figure  5.5  Photographs  of  Dyelines  as  a  Function  of  Depth,  H.  T  =  12.1  s;  6  =  3;  Qo  =  0;  Qf  =  0.52  rad/s;  t/T  =  2  (Cc 
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This  decrease  is  evident  in  column  1  of  Figure  5 .5  as  a  lessening  in  the  stretching  and  wrap¬ 
ping  of  the  dyelines.  Another  consequence  of  decreasing  depth  was  a  corresponding  de¬ 
crease  in  the  spin-up  time.  This  conclusion  is  not  obvious  from  Figure  5.5;  see  §5.7  for 
further  discussion.  Column  2  of  Figure  5 .3g  shows  that  as  a  result  of  the  decreasing  flow 
speed  and  spin-up  time,  the  shallow  fluid  (H  <  1.0  cm)  spun-up  before  3  closed  circulation 
cells  formed.  This  result  was  especially  noticeable  when  H  =  0.5  cm.  In  that  case,  the  fluid 
spun-up  before  the  cyclonic  cells  moved  into  the  tank.  This  decrease  in  spin-up  time,  Tg, 
is  consistent  with  Greenspan  and  Howard’s  (1963)  spin-up  time  scale,  Tg  —  E"1/2Q_1  = 
H/CvQ)1^2 ,  for  fluid  in  a  circular  cylinder  for  which  the  secondary  flow  induced  by  the  Ek- 
man  layer  on  the  bottom  boundary  is  the  mechanism  for  spin-up  (See  §2). 

We  remark  that  the  points  of  concentrated  dye  in  Figure  5.5h  are  the  original  loca¬ 
tions  of  the  dye.  The  striations  in  the  blue  dye  in  column  2  of  Figure  5.5h  are  the  high-fre¬ 
quency  waves  discussed  in  §4.2.  These  waves  became  visible  after  spin-up  when  the  dye¬ 
lines  were  stationary  and  will  be  discussed  in  more  detail  in  §5.6. 

5.4  Effects  of  Varying  fy 

In  the  experiments  discussed  here,  the  fluid  spun-up  from  rest  when  6  =  3,  H  =  5  cm, 
and  Qf  was  varied  from  0.13  rad/s  to  1.82  rad/s.  Figure  5.6a  shows  that  for  the  slowest 
rotation  rate  examined,  the  fluid  behaved  similarly  to  the  shallow  depth  cases  in  that  the  cy¬ 
clonic  eddies  moved  away  from  the  walls  of  the  tank  before  they  formed  closed  circulation 
cells.  Column  2  of  Figure  5.6a  shows  that  the  fluid  spun-up  before  the  center  eddy  closed 
off.  A  doubling  of  this  rotation  rate  (Figure  5.6b)  resulted  in  dyelines  that  were  similar  to 
those  of  the  comparison  case  described  in  §5.1  and  shown  in  Figure  5.6c.  The  dyelines  at 
spin-up  formed  an  essentially  3-cell  pattern  with  two  secondary  anticyclonic  cells  in  the 
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leading  corners  of  the  cylinder.  This  pattern  became  more  complex  when  the  rotation  rate 
was  six  times  the  slowest  considered  here.  After  spin-up  (column  2  of  Figure  5.6d),  the 
dyelines  were  similar  to  those  for  the  case  of  6  =  2.24,  Qq  =  0,  and  Qf  =  0.52  rad/s  described 
in  §5.2  above;  however,  the  evolution  to  this  pattern  was  very  complicated  and  not  like  that 
of  the  6  =  2.24  case.  The  flow  began  similarly  with  a  central  anticyclonic  cell,  which  was 
quickly  separated  from  the  wall  and  deformed  by  the  rapidly  growing  cyclonic  cells  formed 
in  the  leading  corners  of  the  cylinder.  As  the  cyclonic  cells  moved  away  from  the  walls,  anti¬ 
cyclonic  cells  were  left  at  the  end  walls.  These  grew  and  wrapped  into  the  cyclonic  cells. 
At  the  same  time,  new  cyclonic  cells  were  forming  in  the  leading  corners  of  the  cylinder. 
These  secondary  cyclonic  cells  also  got  pulled  into  the  main  cyclonic  cells.  This  process 
continued  as  new  cells  framed  in  the  comers,  grew,  and  got  pulled  into  the  larger  cyclonic 
cells .  Column  2  of  Figure  5.6d  shows  that  fluid  pathlines  after  spin-up  formed  an  essentially 
5-cell  pattern  —  a  center  anticyclonic  cell  flanked  on  each  side  by  a  dominant  cyclonic  cell 
and  a  secondary  anticyclonic  cell.  The  centers  of  the  cyclonic  cells  were  on  the  diagonal 
toward  the  trailing  comers;  the  centers  of  the  secondary  anticyclonic  cells  were  toward  the 
leading  comers.  Figures  5.6e,  f,  and  g  show  that  the  results  for  Qf  =  1.04, 1.30,  and  1.56 
rad/s,  which  are  about  ten  times  the  slowest  rate,  were  very  similar  to  those  just  described. 
Column  1  of  these  figures  shows  that  with  each  increase  in  Qf  the  transient  flow  at  the  onset 
of  spin-up  became  more  turbulent.  This  turbulence  was  not  evident  after  a  few  rotation  peri¬ 
ods  as  the  dyelines  became  smooth.  Column  2  shows  the  fluid  pathlines  at  spin-up.  The 
results  most  closely  resembled  those  of  vHDD,  when  £2f  =  1.82  rad/s  (see  §5.1).  For  the 
large  rotation  rate,  the  dyelines  at  spin  up  formed  an  essentially  3-cell  pattern.  Unlike  the 
results  for  slower  speeds,  and  like  the  result  of  vHDD,  the  center  cell  was  cyclonic  and 
flanked  by  two  anticyclonic  cells.  However,  unlike  the  vHDD  result,  these  cells  were  not 
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circular;  the  center  cell  was  stretched  out  in  an  S-shape  and  the  two  anticyclonic  cells  were 
elongated  on  one  side.  The  evolution  to  this  circulation  was  very  turbulent  and  complicated. 

5-5  Path  of  the  Eddy  Center 

We  tracked  the  path  of  the  center  of  one  of  the  cyclonic  cells  famed  at  the  comer 
for  three  different  Qf’s:  0.26, 0.52,  and  1.04  rad/s  (6  =  3,  H  =  5  cm,  =  0  for  all  cases). 
Figure  5.7  shows  the  measured  paths.  In  all  cases,  the  path  of  the  center  was  a  non-mono¬ 
tonic  spiral  to  its  endpoint  with  the  center  of  the  eddy  ending  up  more  toward  the  center  of 
the  cylinder  than  would  be  expected  if  the  3  cells  equally  divided  the  tank.  The  paths  show 
that  the  eddy  center  was  also  off  the  y -center line  of  the  tank  after  spin-up.  The  centers  of 
the  cells  ended  at  a  location  increasingly  distant  from  the  y -centerline  as  Qf  was  increased. 
This  conclusion  is  not  obviously  supported  by  the  data  for  Qf  =  1.04  rad/s  in  Figure  5.7  be¬ 
cause  the  cell  was  not  tracked  all  the  way  to  spin-up,  but  the  photographs  in  Figure  5.6  clear¬ 
ly  support  the  conclusion.  Figure  5 .8  shows  die  distance  of  the  eddy  center  from  the  comer 
of  the  cylinder  as  a  function  of  time.  As  Qf  increased,  the  time  required  for  the  eddy  center 
to  move  into  the  center  of  the  tank  decreased  from  t/T  ~  5  to  t/T  ~  2. 

5.6  Path  of  a  Separation  Point  of  the  Flow 

We  used  photographs  taken  every  two  seconds  as  described  in  §4.3  to  track  the  path 
of  a  separation  point  of  the  flow.  The  separation  point  formed  along  the  long  wall  (and  stayed 
on  the  wall)  where  the  flow  from  the  cyclonic  eddy  converged  with  the  flow  from  the  anticy¬ 
clonic  eddy,  causing  a  portion  of  the  fluid  to  move  toward  the  wall  where  it  separated  into 
two  small  counter-rotating  eddies  as  discussed  in  §5 . 1 .  Figure  5 .9  shows  the  distance  of  the 
separation  point  along  the  x-endwall  from  the  corner  of  the  cylinder  as  a  function  of  time. 
Once  the  tank  began  spinning,  the  separation  point  moved  smoothly  down  the  side  of  the 


Figure  5.7  Path  of  the  Center  of  the  Cyclonic  Eddy.  H  =  5  cm;  Qq  =  0;  6  =  3;  E_1/2  =  36; 

Qf  =  (a)  0.26,  (b)  0.52;  ( - )  Observer  1,  (-  -)  Observer  2,  (c)  1.04  rad/s.  The  intersection 

of  the  x  =  L/6  and  y  =  W/2  indicates  where  the  center  of  the  eddy  would  be  if  the  3  cells 
divided  the  tank  equally.  (The  l-l  shows  the  error  of  measurements.) 
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Figure  5.8  Distance  from  the  Cylinder’s  Comer  to  the  Center  of  the  Cyclonic  Eddy  as  a 
Function  of  Time.  T  =  12. 1  s;  Qo  =  0;  6  =  3;  E"1/2  =  36;  and  Qf  =  (a)  0.26,  (b)  0.52, 
- Observer  1, - Observer  2,  (c)  1.04  rad/s. 
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Figure  5.9  Path  of  a  Separation  Point  of  the  Flow  along  the  x-endwall.  T  =  12.  Is;  Qq 
=  0;  6  =  3;  E~l/1  =  36;  L/H  =  5.6;  and  Qf  =  0.52  rad/s.  The  dashed  line  indicates  1/3  of 
the  distance  along  the  wall. 

cylinder.  When  it  reached  a  position  slightly  less  than  1/3  of  the  distance  along  the  wall,  it 
oscillated  on  the  wall  about  that  location. 

5.7  Post  Spin-up  Phenomena 

Several  interesting  flow  patterns  were  observed  at  the  latter  stages  of  the  spin-up  oro- 
cess.  We  believe  these  phenomena  occurred  after  the  fluid  achieved  solid  body  rotation  be¬ 
cause  they  occurred  several  revolutions  after  visible  changes  in  the  dyelines  were  noted. 

The  first  phenomenon  observed  was  the  striations  in  the  fluid  dyelines,  wLkh  be¬ 
came  visible  after  spin-up  in  the  shallow  depth  cases  as  seen  in  Figure  5. 10.  These  striations 
were  predominantly  seen  toward  the  ends  of  the  cylinder.  They  are  evidence  of  the  high-fre¬ 
quency  surface  waves  discussed  in  §4.2 


Figure  5.10  Surface  Waves  Evidenced  in  the  Dyelines  after  Spin-up.  T  =  12.1s;  Qq  =  0; 
Qf  =  0.52  rad/s;  b  =  3;  H  =  1.0  cm;  and  t/T  =  6. 


Figure  5.11  Post  Spin-up  Dyeline  Instability.  T  =  12.1s;  Qq  =  0;  Qf  =  0.52  rad/s; 
6  =  2.5;  H  =  5.0  cm;  and  t/T  =  16. 
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The  second  phenomenon  observed  was  the  wavy  pattern  in  the  dyelines  seen  in  Fig¬ 
ure  5. 1 1 .  This  pattern  was  first  discernable  in  the  dyelines  at  the  cylinder’s  center  and  later 
evident  at  the  end  walls.  The  third  phenomenon  noted  was  the  small  vortex  formation  as  seen 
in  Figures  5.12  and  5.13.  These  eddies  formed  first  in  the  center  of  the  cylinder  and  later 
spread  throughout  the  fluid.  Figure  5. 13a  shows  a  close-up  of  the  eddies.  The  small  count¬ 
er-rotating  eddies  behind  the  large  cyclonic  cell  are  reminiscent  of  a  K  arm  an  vortex  street 
formed  by  flow  past  a  cylinder  at  high  Reynolds  numbers.  It  is  likely  that  the  cause  of  the 
second  and  third  phenomenon  is  the  difference  in  surface  tension  and  density  of  the  water 
and  the  dye. 


Figure  5.12  Post  Spin-up  Vortex  Development.  T  =  3.0  s;  Qq  =  0.78  rad/s;  Qf  =  2.07  rad/s; 
6  =  3;  H  =  5.0  cm;  and  t/T  =  (a)  58.25,  (b)  52.25. 


Figure  5.13  Post  Spin-up  Vortex  Development.  T  =  12. 1  s;  Qq  =  0;  Qf  =  0.52  rad/s; 
6  =  3;  H  =  2.0  cm;  and  t/T  =  9.75. 
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Chapter  6 

CONCLUSIONS 

The  presence  of  corners  in  the  rectangular  cylinder  caused  the  formation  of  eddies 
in  the  flow  as  the  fluid  was  accelerated  from  rest  to  solid  body  rotation.  The  number  of  ed¬ 
dies,  as  well  as  eddy  size,  position,  and  rotation  rate,  were  found  to  be  dependent  on  6,  Of, 
and  H.  We  evaluated  our  experiments  against  a  comparison  case  for  which  8  =  3,  Qf  =  0.52 
rad/s,  and  H  =  5  cm.  Horizontal  dyelines  in  this  experiment  formed  a  3-cell,  cyclonic-anti- 
cyclonic-cyclonic  (-  +  -)  pattern.  The  cells  were  near-circular.  Small  secondary  anticy- 
clonic  eddies  were  observed  in  the  trailing  comers  of  the  cylinder. 

As  8  was  decreased  from  3  to  2,  the  size  of  the  secondary  anticyclonic  cell  increased. 
When  8  <  2. 14,  the  flow  evolved  through  3  circular  cells  into  an  arrangement  of  5  non-cir¬ 
cular  anticyclonic  and  cyclonic  cells  that  were  not  aligned  along  a  horizontal  axis. 

As  H  was  decreased  from  5  to  0.5  cm,  a  corresponding  decrease  in  the  spin-up  time 
was  observed.  This  confirmed  the  importance  of  the  Ekman  number,  E = v/(QfH2),  as  a  mea¬ 
sure  of  the  spin-up  time  scale.  A  decrease  in  the  rotation  rate  of  the  cyclonic  eddies  was  also 
observed  as  H  was  decreased. 

As  Qf  was  increased  from  0. 1 3  to  1 .82  rad/s,  the  centers  of  the  cyclonic  cells  became 
more  diagonalized  and  the  tendency  of  the  cyclonic  cells  to  travel  to  the  center  of  the  cylinder 
increased.  These  tendencies  were  confirmed  by  the  measurements  of  the  cyclonic  eddy  paths 
discussed  in  §5.5.  For  a  large  enough  Qf,  the  cyclonic  cells  joined  in  the  center  of  the  cylin¬ 
der,  so  that  the  final  pathline  of  the  flow  was  like  that  of  vHDD.  Comparison  with  the  experi¬ 
ments  of  vHDD  showed  that  the  length  of  the  cylinder,  in  addition  to  Qf,  was  an  important 
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factor  affecting  the  tendency  of  the  cyclonic  eddies  to  move  into  the  center  of  the  tank.  Thus 
El  =  v/(QfL2)  seems  to  be  an  important  scale  that  is  relevant  to  the  evolution  of  the  pathlines 
during  the  spin-up  process. 

Several  interesting  post  spin-up  phenomena,  including  dyeline  instabilities  and  vor¬ 
tex  development,  were  observed.  It  is  likely  that  these  instabilities  are  due  to  the  differences 
in  fluid  properties,  such  as  density  and  surface  tension,  of  the  dye  and  the  water. 
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